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ABSTRACT 

Stability of cylindrical interface between two ideal incompressible fluids, in¬ 
cluding the magnetic field, surface tension and gravitational field is studied in 
linear approximation. We found that helical waves arising both in plasma comet 
tails and on the vertical cylindrical water jet in the air are described by the same 
dispersion equation where the comet tail magnetic field plays the same stabilizing 
role as surface tension for water jet. Hence they represent the same phenomenon 
of Kelvin-Hclmholtz instability. Thus, helical waves in planetary and cometary 
magnetotails as well as in astrophysical jets may be simulated in the laboratory. 

The resonance nature of the instability damping is demonstrated. 

Subject headings: instabilities — comets: general 


1. Introduction 


Stabil ity of p l ane in terface be tween two ideal incompr e ssible fluids has been first con¬ 
sidered by iKelvinl (|187ll) (see also (ILandau &; Lifshitzl Il959l: iMilne-Thomsonl Il960h ). 


Small oscillations in such a fluid are always potential in the first approximation. There¬ 
fore, the velocity potential satisfies the Laplace equation, and for perturbations proportional 
to exp [i{kx — cot)] we arrive at the dispersion equation 
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where p is the density, V is the velocity, g n is the normal (to the interfa ce) co mponent 
of the gravitational acceleration, and a is the surface tension. IRavleighl (118921) was the 
first to consider a cylindrical interface but his stability analysis was restricted by ’’varicose” 
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perturbations (in modern nomenclature, sausage-like ones) also proportional to exp [i{kz — 
cot)}, and hence could not describe kink-modes proportional to exp [i{kz + rncp — cut)], m^O. 


More general analysis was required following a discovery of the geomagnetic tail by .Ness 


(119651) . It has been performed in (IMcKeuziei 1 9701 ) for a cylindrical interface in compressible 
plasma in MHD approximation (the dis persion equation in this case is tra nscendental and 
may be solved only numerically), and in (jErshkovich fc Nusinovlll97lL Il972l ) for incompress¬ 
ible plasma with the magnetic field (like the plasma bulk velocity) parallel to the cylinder 
axis z. The solution of Laplace equation proportional to exp [i(kz + rrup — cut)] describes 
helical waves, like kink modes m = 1,2, etc., along with sausage mode m = 0 (the solution 
is single-valued only if m is an integer, positive or negative). If such solutions are possible 
then helical waves may exist in the nature, identified and observed. 


Indeed, in contrast to the Earth magnetic tail, such oscillations have long been observed 
visually in rectilinear comet plasma tails (type I tails). Bessel (1183 6) was the first to describe 
these wave motions in detail. Alfven (119571) assumed that they are MHD waves. But helical 
waves arising due to interface instability are surface waves propagating in both fluids (as 
a whole). Hence, they cannot be Alfven waves as the Alfven velocity in the comet tails 
(where the plasma is heavy) is usually much less than in the neighboring solar wind. The 
quantitative MHD theory of helical wave origin in type I comet tails due to instability 


in plasma cylind er was suggested in (jErshkovich et.al.l Il972l: lErshkovich &; Chernikovl Il973 


Ershkovichl 11980). The corresponding dispersion equation for the model of plasma cylinder 
with radius R is 
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where indices i and e refer to internal and external plasmas, respectively, and the function 
L m = [I[ykr)K m (kr)]/[I rn (kr)K' m [kr)} is taken at the unperturbed interface r = R. I m 
and K m are modified Bessel functions, a stroke denotes the derivative by argument kr. The 
function L m (kR ) is always negative: 0 > L m > —1 (see Figure 1). 


By using the same standard procedure of linearization described in more details in 
(lErshkovich fe Nusinovlll972l: Ershkovichl 119801 ) we may also include the gravitational field g 
and surface tension a, arriving at the following dispersion equations: for cylindrical interface 
between two liquids 
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and for liquid-gas interface 
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( 4 ) 

The difference in the last term between equations above is due to the fact that a = 0 for the 
gas which is located outside the cylinder (r > R). 


2. Discussion 


With 14 = = 0, er = 0 equations (3) and (4) describe stable ’’surface 

Alfven waves” convected with the fluid bulk velocity V: u/k = V ± [(B 2 — L m Bl) / (An pi — 
AnLmPe )} 1 / 2 . Standard expression Va = B/y/Anp is obtained from the expression above with 
P Pi Pei B Bi B e . 


In the short wavelength limit kR 1, L m —> —1, and with g'_= 0, a = 0 equa tion (3) 
reduces to the dispersion equation for plane interface obtained by IS yrovatsk ii (11953 1 within 
MHD for incompressible plasma. With kR 1 and B = 0 equation (3), naturally, reduces 
to the dispersion equation (1). Finally, with B = 0, g n = 0, V e = 0 equation (4) is the 
solution of the dispersion equation (2.1) in (Yarin 2011) for a verti cal cyl i ndric al water jet in 
air. Strictly speaking, the dispersion equation (2.1) for water jet (1 Yarin 20111) was obtained 
for perturbations proportional to exp [i(kz-ojt)}. Such perturbations are plane waves, which, 
in principle, cannot correspond to helical waves. But it is clear that the structure of the 
dispersion equation for cylindrical jet, being dependent on combinations I m , K m and their 
first derivatives, has the same form for all integer values of m. 


As L m < 0, I’ m /I m > 0 both magnetic field and surface tension terms are always 
positive and hence tend to stabilize the interface. The velocity shear term under the radical 
is negative being responsible for the Kelvin-Helmholtz instability. The term proportional to 
g n results in the flute (Raylcigh-Taylor) instability for heavy fluid over the light one, and 
with pi p e , V = 0 it describes stable gravitational waves on deep water: lo 2 = kg. 


External gravitational field always tends to violate the cylindrical symmetry, because 
the external normal to the cylindrical surface changes its sign (relative to the field) at the 
opposite sides of the cylinder. 


There are, however, several cases when g n is small or even vanishes. This is, of course, 
almost rectilinear (an d cy lindrical) plasma tail of comets, where helical waves are visually 
observed (see e.g. (Ershkovich 1980.) and references therein) as well as astrophysical jets, 
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e.g. (Birkinshawj 1996 ). Formation of plasma comet tails is governed by solar wind, and, 
hence, they are almost in antisolar direction, with rather small aberration angle due to 
orbital motion around the Sun. Another example (although hypothetical) is self-gravitating 
astrophysical jet where the gravitational field is radial and, hence, axially symmetric. 

_ St ability of vertical water jets has been studied in laboratories fo r a long time (see, e.g., 

(|Yarinll2011[ iLeibovich fc Stewardsonlll983l : iGallaire fe Chomazl 120031 1 . Let us now estimate 
the minimal velocity shear V m i n (R ) required for instability of vertical water jet in air. With 
Pi p e , B = 0, V e = 0 (air is at rest), g n — 0, V — Vi, equation (4) yields 


y = V ± L m p e V* + ak 1 ^, 
k \JPi V Im 

whence one obtains the instability criterion in the form 


(5) 


\L m \p e V 2 > ak fk 


( 6 ) 


As functions L rn and I' m /I m depend on kR, the value kR is to be estimated. 


Observations in plasma comet tails show that kR ~ 1 (e.g. ( lErshkov ich 19801 ), Table 
II). This fact is not unexpected as the cylinder radius R (more precisely, the circumference 
2ttR) is the only characteristic length scale of the problem under consideration. Observations 
in water jet are also indicative of kR ~ 1 , and the first attempt to explain this phenomenon 
was performed by Rayleigh ( 18921 ) . It is indeed tempting to explain this fact by means of 
behavior of the instability growth rate, 7 = Imca. But this is not so: the function 7 (kR) has 
no maximum with kR = 1. According to equation (5), 
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Hence, 7 oc k^\L m \ if o — > 0 or the velocity of jet, V, is large enough. The factor 
Xy/\L m (x)\ is shown in Figure 2. It is seen that the dependence of the instability growth 
rate Inna = 7 (kR) is monotonous, and has no maximum with kR ~ 1. 

Fluid parameters are varying in broad ranges while kR ~ 1 remains almost the same. 
It means that kR ~ 1 is a geometric characteristic, a peculiarity of a cylinder. Observation 
mentioned above seems not to be explained within a fram ewor k of infinitesim al amplitudes. 
On the other hand, a finite amplitude theory (lErshkovich fe Chernikov 1973 ) shows that the 
critical wave amplitude 5 C of non-linear stabilization, indeed reaches maximum with kR ~ 1 
as 5 r J R = (Ai) -1//2 , and the function Ai (kR) has here the sharp minimum (Figure 4 in 
(lErshkovich 1980) shows the function A m {kR) for m = 0 and m — 1). 
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According to Laplace formula, an additional surface tension under the cylinder interface 
is cr/R, so that with kR = 1 the coefficient a = 74 dyn cm -1 (for a plane water surface in 
air) is double valued, and using equation (7) we arrive at the conclusion that the instability 
arises if V 2 > (2al' m /I m )/(p e R\L m \). 


Finally, with cr = 74 dyn cm x , (h/p P = 770 for the wate r-air interface, \L±\ & 0.75 and 
I'm / 1m ~ 1-25 for kink mode m = ±1 ( flKruskal fe Tucklll958l) . Figure 1) one obtains 


14 


ms 


-ii 


/ 


20 

7? [cm] 


( 8 ) 


Thus, the minimal initial jet velocity of water vertical free fall required for instability is 
4.5 m s , 2.0 m s , and 1.4 m s~ x for the water jet radius R = 1, 5, and 10 cm, respectively. 
The value V rnin happens to be too high in order to observe helical waves in the cylindrical 
vertical jet from a water tap, but, instead, it is possible to observe there stable surface waves 
(i.e. normal modes of oscillations) traveling along the jet both upstream and downstream 
(in the frame of reference moving with the velocity V) when the expression under the radical 
in the equation (5) is positive. This phenomenon is just the same as ripples created by wind 
on the lake surface. 


Flelical waves in plasma comet tails may become visible only when their amplitude 
becomes large enough. As a result, this phenomenon is observed relatively seldom. The 
stability conditions seem to be marginal. As helical waves in comet tail and in vertical water 
jet obey almost the same dispersion equation this astrophysical phenomenon may be (and, 
in our opinion, is to be) simulated in laboratory. 


We did not consider here the effect of finite width, d, of a transition layer between two 
fluids. It is known to be small with kd 1 , and as for cylindrical jet kR ~ 1, this effect 
is negligible if d <C R which is the case for water jet in air and seems to be observed (as 
sharp decrease of brightness) for plasma comet tails. Model of cylindrical comet tail with 
transition layer of finite thickness d is treated by Chen & Liu (119821 ). 


3. Resonance damping of helical waves 


Landau! ([1.944) found a sharp decrease of the Kelvin-Helmnoltz instability growth rate 


when the phase velocity of surface wave, Reu/k is approaching the acoustic velocit y, c, w ith 
full d amping 7 = Imcn = 0 when the phase velocity r eaches \[2c (see also (Eiandau &; Lifshitz 
19591 ). ch.9, § 84). A similar effect was described in (Ray & Ershkovich 1983) for magnetoa¬ 


coustic velocity. We believe that there is a simple explanation: when the phase velocity of 
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surface wave, Rec o/k, approaches the characteristic velocity of normal mode of oscillations 
in the fluid, a resonance arises, and stable hydrodynamical or MHD modes are generated in 
the whole fluid volume. But the Kelvin-Hclmholtz instability of the tangential discontinu¬ 
ity is a surface phenomenon, with the amplitude of perturbation decreasing (in the plane 
exponentially) away from the interface. Thus, these waves, in some meaning, are 


case 


two-dimensional, with relatively restricted stock of kinetic energy, supplied by the velocity 
shear. When this energy is transferred from 2D to 3D space generating stable waves ev¬ 
erywhere, the energy stock is rapidly exhausted, and instability is damping. But if so, the 
same phenomenon should exist in incompressib le plasma while reaching the Alfvenveloc ity, 
oo/k ~ Va- Indeed, this effect was described in (IRav fe Ershkovichlll983t lLau fe Liulll98ll) for 
plane interface. 


Let us consider a cylindrical plasma jet (with the velocity V) immersed into plasma at 
rest, with the same parameters, i.e. Pi = p e = p, B t = B e = B , and V e — 0 (alternatively, 
we may choose the frame of reference where V e = 0). Then equation (2) yields the phase 
velocity Reoo/k = V/{1 — L m ), where V = Vi is the velocity jump. The radical in equation 
(2) vanishes, and the interface becomes stable with V = (1 — L m )V^/|L m |, and the phase 
velocity ou/k = V/{1 — L m ) = VA/\L m \ = 1.15Va for kink mode m = 1 ,L m = —0.75. For 
plane interface L m = — 1 , hence V = 2 Va, and ou/k = Va- 


A possibility of resonance generation of Alfven waves in the whole volume of fluid 
with oo —y kVA seems to be obvious. The fact that the magnetized shear layer is stable 
i f its Alfven speed is greater than half th e velocity jump across the interface was found in 
(IRav fe Ershkovichlll983nLau fe Liull 198 lh (unfortunately, the resonance nature of the i nsta- 
bility damping has not been mentioned therein). According to IRav & Ershkovi ch (11983 ). the 
interface remains stable (despite the growing flow velocity V) when V >2c and Va > c. As 
the phase velocity oo/k = V/2 , these conditions may be rewritten in the form oo/k > c <Va 
, which agrees with the resonance scenario above: the instability ceases because the energy 
supplied by the velocity shear transfers (due to the resonance) for excitation of normal modes 
of the fluid oscillation, first, of sound waves (as c < Va), and then, of MHD waves (Alfven 
and magnetosonic). 


Similar resonance damping occurs with unstable capillary waves. In order to demon¬ 
strate this effect, let us assume that the liquid in a cylindrical jet (moving with the speed 
V) has almost the same density, p tl as the ambient liquid at rest, i.e. \pi — p e \ -C p^ e = P- 
Then equation (3), with B = 0, g n = 0 yields 
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The radical vanishes if 


and the phase velocity 
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Hence, for a kink mode m = 1, L m = —0.75 one obtains cu/fc = 0.87 \Jak/p. For the plane 
interface L rn = — 1 we find u/k = \Jak/{2p) (which, naturally, may be obtained directly 


from equation(l)). The classi cal value for stable capillary waves is c o/k = \J akj p (Ke lvin 


1871; Landau & Lifshita 1959). 


A small region of the cylindrical interface may be considered as plane for perturbations 
with kR 1, and for plane case the dispersion equation for perturbations ~ exp [i(kr-wt)] 
depends on scalar products kV and kB. This means that there are always directions along 
which the stabilizing role of the magnetic field becomes negligible. As short wavelength 
perturbations may propagate in all directions a tangential discontinuity always remains un¬ 
stable. But this is not the case for helical waves propagating along the cylinder axis. 


4. Conclusion 


Dispersion equations (3) and (4) describe rather broad class of hydrodynamical and 
MHD instabilities and normal modes of oscillations of the cylindrical interface between two 
fluids, started with Alfven waves and gravitational waves on deep water to flute and Kelvin- 
Hclmholtz instabilities in planetary and comet tails and water jets in air (including also 
capillar instability in liquids). Although they were obtained in linear approximation (and 
hence each of these effects may be studied independently) the stability criterion Imw = 0 de¬ 
pends on the balance of all the relevant terms under the radical. This balance is particularly 
important under marginal stability conditions when only their sum is indicative of stability 
or instability of the interface. 


We also drew attention to the fact that the instability growth rate obtained in linear 
approximation cannot explain the preferential generation of modes with kR ~ 1. In partic¬ 
ular, this fact refers to helical wa ves observed visually in comet plasma tails . At t he same 
time, finite amplitude treatment flErshkovich fe Chernikovl Il973l: lErshkovichl Il980l ) seemed 
to explain these observations. 


Both Alfven and capillary waves arising due to Kelvin-Helmholtz instability on the 
cylindrical interface have been considered. We arrived at the conclusion that sharp damping 






















of these helical waves occurs when their phase velocity approaches the characteristic velocity 
of normal modes of oscillation, so that it has resonance nature. 

Finally, we found that helical waves both in plasma comet tails and in vertical cylindrical 
water jet in the air are governed by almost the same dispersion equation (which means that, 
in fact, we deal with the same phenomenon). This fact allows us to suggest an idea of 
laboratory simulation of helical wave generation in cometary and planetary magnetotails as 
well as in astrophysical jets by using vertical water (or any other suitable liquid) jet. 
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Fig. 1.— T he function L rn (kR) wit h m = 0 and m 
according to Ershkov ich & Chern ikov (119731 )). 
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Fig. 2.— The function xy/\L m (x)\ with m = 0 (solid line) and m = 1 (dashed line), x = kR. 








